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Abstract: We consider a model of preheating where the coupling of the inflaton to the 
preheat field is modulated by an additional scalar field which is light during inflation. We 
establish that such a model produces the observed curvature perturbation analogously to 
the modulated reheating scenario. The contribution of modulated preheating to the power 
spectrum and to non-Gaussianity can however be significantly larger compared to modu- 
lated perturbative reheating. We also consider the implications of the current constraints 
on isocurvature perturbations in case where the modulating field is responsible for cold 
dark matter. We find that existing bounds on CDM isocurvature perturbations imply that 
modulated preheating is unlikely to give a dominant contribution to the curvature pertur- 
bation and that the same bounds suggest important constraints on non-Gaussianity and 
the amount of primordial gravitational waves. 
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1 Introduction 

One of the key issues for inflationary models is the question of how the inflaton energy is 
converted to radiation of ordinary particles. In the conventional theory of reheating [1-3] 
the inflaton decay is perturbative with an effective decay width T. In a typical model T 
is much smaller than the Hubble parameter H at the end of inflation so that at first the 
inflaton loses more energy to the expansion of the universe than is transferred to its decay 
products. Once the Hubble parameter decreases to H = T, reheating becomes efficient 
and the inflaton energy is transferred to relativistic particles. In the so-called modulated 
reheating scenario [4-6] the inflaton decays into radiation also in the perturbative regime 
but with a width that is dependent on another, spatially fluctuating scalar field. This 
modulating field is assumed to achieve perturbations during inflation, giving thus rise to a 
space-dependent decay rate. 

In models with several scalar fields the decay of the inflaton can also take place in a non- 
perturbative regime through parametric resonance [7-9] (for recent reviews see [10, 11]). 
In this case, the inflaton energy is converted to particles in violent bursts, followed by a 
later thermalization of the decay products. 

If the coupling of the inflaton to the preheat field depends on yet another light field, 
we would have modulated preheating: inflaton would decay non-perturbatively but the 
parametric resonance would be space-dependent. It has been argued that the modulation 
of effective, field-dependent coupling constants follows naturally from string theory, and the 
implications for preheating were discussed in [12] and [13]. Curvature perturbations from 
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modulated preheating were considered in [14, 15] but without the effects of backreaction. 
Modulated perturbations in the context of preheating in hybrid inflation were studied in 
[16], and in the context of trapped inflation in [17]. 

If several fields are dynamical during inflation, in addition to the adiabatic perturbation 
there will also arise an isocurvature perturbation. These can feed the curvature perturba- 
tion [18-20], but even if the additional fields do not evolve during inflation, the isocurvature 
perturbations from those fields can be converted into the curvature perturbation once in- 
flation is over. This is the case for the curvaton scenario [21-24], where the curvaton field 
may come to dominate the energy density of the universe after inflation (see [25] for a 
recent review) . This also holds for the modulated reheating scenario [4-6] where the decay 
width of the inflaton depends on an additional field. In addition to isocurvature pertur- 
bations, another signature of multifield scenarios could be non-negligible non-Gaussianity 
that might be observable already by the Planck satellite. 

In the present paper we consider an extension of the idea of modulated reheating to 
the theory of parametric resonance. We study the case where the coupling of the inflaton, 
denoted as <j>, to the preheat field \ is modulated by its dependence on an additional scalar 
field, denoted as a. One possibility is that the modulating field a could be dark matter, 
and we investigate the implications of isocurvature constraints on such a model. We also 
study the non-Gaussianity and the dependence of the non-linearity parameter /jvl on the 
modulating field. 

The outline of the paper is as follows. In Section 2 we briefly review the modulated 
reheating scenario and the theory of parametric resonance in an expanding universe. In 
Section 3 we consider the effects of the modulated coupling in preheating and compute 
the power spectrum and non-Gaussianity. In Section 4 we consider the constraints on the 
model from observational bounds on CDM isocurvature perturbations if the modulating 
field is responsible for dark matter. We conclude in Section 5 with a discussion of the 
results. 

2 Background 

Let us begin by setting the stage for modulated preheating by introducing the notation and 
reviewing the ingredients needed for discussing modulated decay and parametric resonance. 

2.1 Modulated reheating 

In the modulated reheating scenario, the decay width is inhomogeneous due to its de- 
pendence on a field a which was light during inflation and therefore acquired an almost 
scale-invariant spectrum of perturbations. Between the end of inflation and the time when 
H = r(u) the inflaton oscillates in a harmonic potential and the evolution of the universe 
is matter-like, H oc a -3 / 2 , where a is the local scale factor. After the energy is trans- 
ferred from the inflaton to other particles at H = T, the universe becomes locally radiation 
dominated, H oc a~ 2 . Due to the inhomogeneity of the transition some parts of the uni- 
verse will have expanded more than others resulting in density, or alternatively, curvature 
perturbations on the hypersurface of constant density. 
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As is well known, the curvature perturbation on uniform density hypersurface can be 
denned in a gauge-invariant manner as [26, 27] 

C = ~^^-X (2-1) 
3(1 + w) p 

Similarly, for any non-interacting component i we can define the curvature perturbation 
on the hypersurface where the energy density of that component vanishes as 

^ = -^ + 7777^— , (2-2) 
3(1 + Wi) pi 

so that the total curvature perturbation is given by 

t-r 1 {l + w)p 

For practical purposes, in modulated reheating we may treat the different regions as inde- 
pendent FRW- universes and make use of the 5N formalism [19, 28], where the curvature 
perturbation £ is given by the difference of the expansion rate of separate universes with 
£ = 5N(a). The local number of e-foldings after the end of inflation is given by 

( a-end ) 3 ( H en d ) 2 ( T ) ^ ^ 

For modulated reheating the curvature perturbation on uniform energy hypersurfaces is 
then 

IF' g' 

Cmr = SN = N a Sa* = --— £cr* = -— £cr*, (2.5) 

6 T 3g 

where 5a* refers to field perturbation at the time of horizon crossing, and the last equality 
comes from assuming Tag 2 with g = g(cr) a coupling constant. Here and in what follows 
is shorthand for dN/dcj) for any field (p. 

2.2 Parametric resonance in expanding universe 

The simplest example of an effective potential giving rise to a parametric resonance is 

V = l -m 2 4> 2 + \g 2 4>\ 2 , (2.6) 

where <j) is the inflaton and x 1S t ne preheat field. The theory of parametric resonance in 
an expanding universe for such a quadratic model has been studied in detail in [9]. After 
the end of inflation, the inflaton oscillates in a harmonic potential around the minimum, 
and the evolution of the universe is matter-like. Every time the inflaton goes through the 
minimum of the potential the effective mass of x becomes zero, and x particles are copiously 
produced in a resonant regime, which is characterized by its resonance parameter 

— M 

where $o is the amplitude of the inflaton at the end of inflation. 



- 3 - 



The number density of x particles during the resonance grows as (see the Appendix): 



n x (t) - ^L e^, (2.8) 

where /U is the effective Floquet exponent describing the growth rate of the fluctuations. 
The growing number of x particles changes the masses of the fields, causing a backreaction 
which eventually shuts down the resonance. Backreaction becomes important when n x (t) ~ 
m 2 Q(t)/g [9]. Comparing with eq. (2.8) then gives the time at which this happens as 



mil = t~ In 
4/i 



A 



2 m f 3Hq^ (mil) /x 



$n I 2m / g 5 



(2.9) 



where we 1 find A = 16\/2tt 3 . Hq is the Hubble parameter at the beginning of the oscillation 
period so that a(t) = (1 + 3H t/2) 2 / 3 ~ (3# i/2) 2/3 at late times. This first stage of 
preheating is followed by a second stage where the effects of backreaction must be taken 
into account; however, as discovered by Kofman et al. [9], the second stage lasts a short 
while compared to the expansion rate of the universe. Hence we can consider the time ii 
to mark the end of preheating. 

The produced particles are not in equilibrium. Lattice simulations show however that 
the effective equation of state w = (p) / (p), where p and p are the pressure and energy 
density respectively, rises rapidly after a matter-like phase and levels off around a value 
that is typically somewhat below the radiation- like value w = 1/3 [13, 29-31]. 



3 Modulated preheating 
3.1 Curvature perturbation 

If the coupling g of the inflaton to the preheat field x m e Q- (2-6) depends on an additional 
scalar field a which is light during inflation, g = g(cr), we expect that entropy perturbations 
in this field will cause a spatial variation in the process of preheating in analogy to the 
modulated reheating scenario. In the modulated reheating scenario, curvature perturba- 
tions are generated because of inhomogeneities in the time of transition of the equation of 
state from matter- like (w = 0) to radiation-like (w = 1/3). More generally, the scale factor 
grows as a oc £ 2 / 3 ( 1 + u ') so that if we have an instantaneous transition at some time ii from 
w = Wi to w = Wf, the number of e- foldings is 



/ai\ a\ 2 l 1 \ ii 2 1 t , s 

N = In — + In — = In — + In - (3.1 

\ai J \ai ) ?>\\+Wi 1 + WfJ t 31 + u>/ t 

If the transition is inhomogeneous, that is, if ii = ti(cr), curvature perturbations will be 
generated due to different expansion rates in different parts of the universe: 



1 Using the equations of [9] gives A — 64-7T 3 . However, it seems to us that [9] has an error in their steepest 
descent formula. We have also used a different estimate for the location of the maximum in the Floquet 
index. See Appendix A for details. 
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Cmp = SN 



2 Wf dlnti g' 



5<J: 



(3.2) 



31 + Wf ding g 



where Wi is assumed to vanish since the evolution during the period of oscillation of the 
inflaton is matter-like. Here ' denotes the derivative with respect to the field a, and the 
subscript "MP" signals that the curvature perturbation contribution is due to modulation 
of preheating. In fact, as found by Podolsky et al. [13], the final value of w as well as 
the shape of the transition also depend on g and therefore the last term in eq. (3.1) would 
also contribute to the curvature perturbation. Detailed investigation of this term would 
require a comprehensive study of thermalization that can only be done through extensive 
lattice simulations. Thermalization after preheating in a model with a quartic potential 
has been studied in [32, 33]. Here we assume that this effect is subdominant and that the 
subsequent evolution is independent of the coupling. 

If the process of preheating terminates before backreaction becomes important, t\ oc g 
[9] , and the curvature perturbation is of the same order (but opposite sign) as that produced 
in the modulated reheating case. If backreaction shuts down the resonance then the end 
of preheating is given by eq. (2.9) and we get: 



The duration of preheating, mt\, is typically of the order of 10 2 and typically ft ~ 0.13. 
If the term din fi/dln g is negligible, then the amplitude of the curvature perturbation 
is | Cmp | < 10 _1 |Cmr|> where Cmr is the curvature perturbation in modulated reheating 
one would obtain with the same set of parameters. As a consequence, the contribution of 
modulated preheating to the power spectrum in this case is less than 10 -2 relative to the 
corresponding modulated reheating scenario. 

However, we find by solving numerically the equation of motion for the preheat field 
that the typical value of din [i/ ding is O(10). Thus this term can be dominant. In this 
case | Cmp | ^ 5 |Cmr| and the contribution of modulated preheating to the power spectrum 
can be significantly larger than in the modulated reheating scenario with a corresponding 
coupling profile. In Figure 1 we compare the curvature perturbation produced by modu- 
lated preheating to the curvature perturbation produced by the corresponding modulated 
reheating model in the broad resonance regime, go ^ 1> where the interesting effects can 
be found. We take qo ~ 10 4 , and as can be seen in Figure 1, for a range in the resonance 
parameter qo of this order modulated preheating leads to a relatively high contribution. 

3.2 Non-Gaussianity 

The primordial curvature perturbation is given by the combined contributions from infla- 
tion and modulated preheating, 



3 lnii 
d In g 



4/imii — 1 5 dlnfi 



A/imti — 3 [4/imti — 1 ding 



(3.3) 



C = NfiSfa + N a 5a, 



(3.4) 



and the power spectrum, defined by (k)5(k — k') 




(CkCk')> is g iven b Y 
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Figure 1. Comparison of the curvature perturbation from modulated preheating with curvature 
perturbation from modulated reheating as a function of the resonance parameter q a , for Wf = 
0.25, Hq = 0.25m. We have estimated d\n[i/d\ng by first solving the time at which the number 
of produced particles exceeds the limit where backreaction becomes important for each go (to the 
accuracy of half an inflaton oscillation) and then solving the derivatives of [i at that time from 
eq. (2.8). 



r C (k) = (N} + N*) (—J . (3.5) 

Defining parameters 



A = -iV CT cr* = — — , £ = 8e* — M , (3.6) 

2 l + Wf 3g ding \ a * J 

we can write = (1 + ^A 2 )Pi n f with V m i = N^(H/2tt) 2 being the contribution of the infla- 
ton to the power spectrum. Here we have adopted the generic prediction N 2 = l/2e*Mpj 
where e* is the slow-roll parameter describing the deviation from deSitter space evaluated 
at the time of horizon crossing. We assume that the dependence of the coupling on the 
field a is of the form 

g 2 (a)=g 2 (l + a 2 /M 2 ) , (3.7) 

where M is some energy scale. According to our numerical estimates for d In \xjd In g, the 
parameter A is at most 0(1) and so modulated preheating can dominate the curvature 
perturbations only if £ is large, that is, if the amplitude of the field a during inflation is 
sufficiently below Planck scale. 

Non-Gaussianity is given by the expression 



■Jnl = r at Ar c i2 ' ( 3 - 8 ) 



6 = N a N b N ab 
5 lNL [N C N C 

where we have ignored the loop corrections. The part coming from the inflaton is suppressed 
by slow-roll parameters and the only potentially large contribution is N a N a N aa / (N a N a + 
N^N^) 2 . In our parameterization N 2 /N? = £A 2 and 
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N a =g' 



,8N 
dg ' 



N„ 



,8N 



+ g' 



r2 



d 2 N 
dg 2 



(3.9) 



so that we can express non-Gassianity as 



MP 
NL 



i e 2 A 4 



2(1+£A 2 ) 2 



g"<?* 



g'a* 



n 


i 








(T) 



(3.10) 



If a/M <C 1 the first term in the square brackets should dominate. In fact this term 
is of the same form as that obtained from the modulated reheating scenario. Therefore we 
can relate the non-Gaussianity produced by these two processes in this limit: 



MP 
NL 



2w 



f 



' 'dhi^V MR 

JNL ■ 



(3.11) 



1 + Wf ) \d\ag 

This relation of plotted in Figure 2 for different values of the resonance parameter with the 
same estimation for din /i/dlng as before. 
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Figure 2. Comparison of non-Gaussianity from modulated preheating and non-Gaussianity from 
modulated reheating for Wf = 0.25, H = 0.25m in the limit a/M <C 1. Here the same estimation 
for dlnfi/dhig is used as in the previous figure. 



In the oppsite limit, a/M 3> 1, we can can also relate the two non-Gaussianities if £ is 
small by solving £ in terms of fjf^- We find 

w f W dm/A 3 / w f V / dln/A 2 d 2 ln/i j MR 
1 + Wf) \d\ng) \l + Wf) \dlng J d(\ng) 2 J lNL ' 

and the result is plotted in Figure 3. If £ is large, on the other hand, the fj^jj ~ 5/2 in the 
limit a/M 3> 1 whereas non-Gaussianity from modulated preheating can still be large as 
long as £A 2 > 1. It is clear that in all cases modulated preheating can produce much more 
significant curvature perturbations and non-Gaussianities than the modulated reheating 
with the same coupling profile. 
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Figure 3. Comparison of non-Gaussianity from modulated preheating and non-Gaussianity from 
modulated reheating for wj = 0.25, H = 0.25m in the limit a/M S> 1. dln/j,/d\ng and 
d 2 ln/i/<9(ln g) 2 were again estimated by solving them from the equations of motion after 15 in- 
flaton oscillations. 



4 Isocurvature perturbations 

Let us now assume that the modulating field a is responsible for cold dark matter (CDM) 
whereas the preheat field x wm eventually give rise to radiation that consists of the Stan- 
dard Model degrees of freedom. 

If a is subdominant at the beginning of the radiation domination era when both <fi 
and x have decayed, we should have £ ~ £ r , where r refers to radiation. The curvature 
perturbation in CDM is given by Qcdm = Co-- On a hypersurface of uniform radiation 
density the isocurvature perturbation between the CDM field a and radiation is from 
eq. (2.2) 

5 CT = 3(C CT -Cr) = T ^^- (4.1) 

1+ W a P(j 

Assuming the potential for a is harmonic with V(a) = ^m^a 2 and that the mass <C H 
throughout inflation and reheating, a will start evolving only once the Hubble parameter 
has decreased to a value comparable to the mass. If a is the origin of CDM, we can write 
the CDM isocurvature perturbation as 

Scdm — 25a* /a* , (4.2) 

where 5a* is the field perturbation on the initial flat hypersurface, as usual. 

Using the parametrization presented in the previous section in (3.6) we can write the 
power spectra as 

V ( = (1 + £A 2 )Pinf, P s = CPud, C cs = -(XV ini (4.3) 

where and Vs are the curvature and isocurvature spectra respectively, and C^s the 
correlation spectrum. In the above we assume that the perturbations in different fields are 
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uncorrelated at the horizon exit. The isocurvature fraction a and the correlation angle A 
are given by the expressions 



a = 1 s - 5 cos 2 A = - ^ u 4) 

Pc + Ps ~ l + e(l + A 2 )' + 1 ] 

The parameter a has been constrained from the CMB and large scale structure obser- 
vations to be small [34-48]. If the contribution from modulated preheating dominates, i.e. 
£A 2 S> 1, we have a = 1/(1 + A 2 ) and the perturbations are completely (anti-)correlated. 
For completely correlated perturbations 2 a < 0.011 (0.0047) from WMAP7 (WMAP7 + 
BAO + SN) data at 95% confidence limit [46, 47]. This requires A > 10 (15) which is not 
achievable according to our numerical estimations for d In \ijd In g. Therefore the curvature 
perturbation cannot be solely due to this effect. The constraint on anti-correlated cased is 
of the same order as on the correlated case [48]. 

If modulated preheating does not significantly contribute to the curvature perturbation, 
£A <C 1, perturbations are uncorrelated and we have a ~ £/(l + £).The constraints on 
isocurvature for the uncorrelated case are less stringent: a < 0.13 (0.077) [46, 47], and we 
obtain the constraint £ < 0.15 (0.084). 

The bounds on isocurvature perturbations imply limits on the amount of non-Gaussi- 
anity that can be produced. In Figure 4 we show the regions of the parameter space giving 
significant non-Gaussianity along with the constraints from isocurvature 
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Figure 4. Non-Gaussianity and isocurvature bounds for qo — 10250. The blue, the violet and the 
brown regions correspond to Jnl > 5,25, and 50 respectively. The dashed black line represents 
Jnl — 100. The red curves show the isocurvature constraints with the solid, dashed and dotted 
curves giving respectively a = 0.13,0.077 and 0.01. 



2 Note that our sign convention differs from that of [46, 47] and our correlated case corresponds to their 
anti-correlated case and vice versa. 
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and in Figure 5 the dependence between non-Gaussianity and the amount of isocur- 
vature for go = 10250, which is the value for which preheating lasts about 15 inflaton 
oscillations with $> /m = 10 6 . 
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Figure 5. The region of allowed non-Gaussianity and isocurvature for qo — 10250. 

The general trend is the same for other values. The result is similar to that of Takahashi 
et al. [49] who investigated modulated reheating with gravitino dark matter. They found 
that significant non-Gaussianity requires large isocurvature perturbations. 

We note that if the value of the modulating field a at horizon crossing is required to 
be below the Planck mass, then from the definition (3.6) we obtain the condition r* < 2£ 
where r* is the tensor-to-scalar ratio at horizon crossing. Thus 2£ = r max gives an upper 
bound for primordial tensor perturbations. In particular, if we fix /jvx, we can obtain a 
relation between r max and a. This relation is displayed in Figure 6. Detection of non- 
Gaussianty implies bounds on the amount of isocurvarture and gravitational waves. For 
large values of A the constraint on gravitational waves becomes even more stringent because 
r = r*/(l + r*A 2 /2). 

5 Conclusions 

Motivated by the modulated reheating scenario we have investigated a model where the 
decay of the inflaton is not perturbative but happens because of resonant production 
of particles of the preheat field x which is eventually shut down by backreaction of the 
produced particles. We have studied a simple potential of the type g 2 4> 2 X 2 an d have 
limited our analysis to preheating and have not considered the subsequent thermalisation. 
The spatial modulation of the inflaton coupling g to the preheat field is assumed to be 
due to an additional scalar modulating field a, which is light during inflation. Therefore, 
during inflation the modulating field achieves entropy perturbations which are converted 
to curvature perturbations at preheating. We find that the contribution due to modulation 
to the resulting power spectrum and non-Gaussianity can be much more significant than 
if one assumes modulated reheating with the same set of parameters, as was discussed in 
section 3.2. 
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Figure 6. The allowed region for r and a for q$ = 10250 with fixed /jvl = 5, 25, 50. Bound on 
isocurvature imply a constraint on primordial gravitational waves. The line terminates for small 
a because small isocurvature implies small non-Gaussianity (see Fig. 5) so that at some point the 
fixed Jnl can no longer be achieved. 

We also considered the possibility that the modulating field is CDM while the radiation 
in the universe is due to the dynamics of the preheat field. In this case there will be an 
isocurvarture perturbation in CDM which is due to the perturbation in the modulating field 
a. We find that current observational constraints on isocurvature imply that the primordial 
curvature perturbation cannot be solely due to modulated preheating. We also find that 
isocurvature bounds can put severe limits limits on the amount of non-Gaussianity and 
primordial gravitational waves that can be observed. We have not performed a complete 
scan through parameter space but rather presented in Figure 5 an example with the reso- 
nant parameter qo = 10250 that demonstrates the general trend in modulated preheating, 
the presence of non-Gaussianity implies a constraint on the amount of isocurvature. 
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A Estimating the behavior of n 

Here we estimate the magnitude of the terms din /i/dlng and d 2 \nfi/d(lng) 2 from the 
preheating dynamics. The inflaton evolves as 4>(t) = sin(mt), where the amplitude 
decreases as <£(i) = <3?o a ^ 2 (^)- The equation of motion for the x held (in Fourier space) 



is then 




X k + u 2 k (t)X k = 0, 
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where we have defined the frequency = k 2 /(m 2 a 2 ) + 4gsin 2 i with the resonance 

parameter q = g 2 ^ 2 /(4m 2 ) . Time is now measured in units of m _1 . The occupation 
number for particles of mode k is given by 



2 



l*k| 



LOT 



+ \x k \ 



(A.3) 



As the innaton goes through the minimum of the potential, the particle density jumps 
rapidly either up or down, depending on the phase acquired since the last oscillation. 
Density is twice as likely to increase than to decrease so that overall grows. The 
behavior of is plotted in Figure 7. 




Figure 7. n k (t) for q = (32tt) 2 , k 2 /m 2 ~ 40 



We define an effective Floquet exponent \ik as 

fi k (t) = T^- t H2n k ). (A.4) 
The number density of x particles is then obtained by integrating over all the modes 

n x (t) = -L__ J tf knk{t) = -±- 3 J dk k 2 e 2 ^. (A.5) 

For late times the integral can be evaluated using the saddle-point (Lagrange) method, 
which gives 



1 

n x (t) 



47r 2 a 3 



mt 



" k 2 e 2^ k mt 



(A.6) 

fc— krn. 



where ' denotes a derivative with respect to k and k m is the value of k that maximizes 
Hk, and fJ-k=k m is the effective Floquet exponent /i used in Section 3. We can estimate 
fx ~ (Ak) 2 where Ak is the width of the resonance band, which can be estimated as 
Ak ~ x/^/^Qo^ 771 (see [9]). Taking the k m to be in the middle of the resonance band, 
k m ~ | y / 2/7T(?Q m, we obtain equation (2.8). 

Figure 8 shows the behavior of /ifc(i) as a function of k 2 /m 2 for qo = (32n) 2 after 15 
innaton oscillations. As qo is changed the positions and magnitudes of the peaks change 
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Figure 8. (if. as a function of k 2 /m 2 for q = (327r) 2 after 15 inflaton oscillations. 



continuously and therefore so does the global maximum. However the change is not mono- 
tonic because as the peak giving the global maximum goes down another peak may grow 
past it and become the new global maximum. Figures 9 and 10 show how the global 
maximum changes with go- 



(*„,/m) 2 




104d0 



Figure 9. The location of the global maximum of /Xfe as a function of go- The dashed vertical lines 
correspond to places where different peaks take over and the places where the maximum is at the 
origin. In the vicinity of these points the analytic estimates are not reliable. 



The analytic estimates given above are not reliable in the vicinity of the transition 
from one peak to the next because at that point there are two separate peaks giving the 
same maximum and the saddle point approximation is not applicable. The estimates are 
also unreliable for the cases where the global maximum is at k = as can be seen from 
eq. (A. 6). For this reason, we numerically obtain the number density n x by calculating 
the integral (A. 5) and use eq. (2.8) as the definition of [i. This way we obtain the correct 
behavior for all values. The behavior of d In \xjd In g and d 2 In fi/ d(ln g) 2 after 15 oscillations 
of the inflaton is shown in Figures 11 and 12. 
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Figure 10. as a function of qo after 15 infiaton oscillations. The dashed red curve corresponds 
to calculating the global maximum of jik and the solid black curve is obtained by using eq. (2.8) as 
the the definition of /i and solving for it from the integrated n x . 



din ji 




Figure 11. d\nn/d\ng calculated after 15 oscillation of the infiaton for q ~ 10 4 . 



cr In [i 




Figure 12. d 2 ln/i/9(ln g) 2 calculated after 15 oscillation of the infiaton for q ~ 10 4 
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